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ABSTRACT 

We present a seismic study of the (3 Cephei star 6 Ophiuchi. Our analysis is based 
on the observation of one radial mode, one rotationally split £ ~ I triplet and three 
components of a rotationally split £ = 2 quintuplet for which the m-values were well 
identified by spectroscopy. We identify the radial mode as fundamental, the triplet as 
Pi and the quintuplet as gi . Our NLTE abundance analysis results in a metallicity and 
CNO abundances in full agreement with the most recent updated sola r valu es. With 
X e [0.71,0.7211] and Z e [0.009,0.015], and using the Asplu nd et al. (|2005l) mixture 



but with a Ne abundance about 0.3 dex larger (Cunha et al. 20061, the matching of 
the three independent modes, enables us to deduce constrained ranges for the mass 
{M = 8.2±0.3 Mq) and central hydrogen abundance {X^ = 0.38±0.02) of 9 Oph and 
to prove the occurrence of core overshooting (aov — 0.44±0.07). We also derive an 
equatorial rotation velocity of 29±7 km s~^. Moreover, we show that the observed 
non-equidistance of the £ = I triplet can be reproduced by second order effects of 
rotation. Finally, we show that the observed rotational splitting of two modes cannot 
rule out a rigid rotation model. 

Key words: stars: early-type - stars: individual: 9 Oph - stars: oscillations - stars: 
abundances - stars: interiors 



1 INTRODUCTION 

For a few years many breakthroughs have been achieved in 
the field of asteroseismology of /3 Cephei stars thanks to 
several observational photometric and spectroscopic multi- 
site campaigns and long-term monitoring dedicated to this 
kind of pulsating B-type stars. The aim is to improve our 
knowledge of their internal structure and more precisely of 
two mixing processes that affect their evolutionary path: 
core overshooting and internal rotation. 

The first detailed asteroseismic modelling was per- 
formed for V836 Cen, which led to constraints on global 
stellar parameters but also on the core overshooting param- 
eter. Moreover the non-rigid rotation of the star was proved 



(Aerts et al. 2003 Dupret et al. 20041. Similar results were 
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afterwards obtained for u Eri (Pamyatnykh et al. 2004 Aus- 
seloos et al. 2004|. Recently, Aerts et al. (2006) gave con- 



straints on the physical parameters of S Ceti thanks to the 
discovery of low-amplitude modes by the satellite MOST. Fi- 
nally, the seismic interpretation by Mazumdar et al. (20061 



showed the occurrence of core overshooting for /3 CMa. 
The derived overshooting parameter values are 0.10±0.05, 
0.05±0.05, 0.20±0.05 and 0.20±0.05 for V836 Cen, u Eri, 
S Ceti and /? CMa, respectively. 

The P Cephei star 6 Ophiuchi was also the subject of 
intensive photometric and spectroscopic observations as de- 
scribed in Handler et al. (20051 and Briquet et al. (2005| 



(hereafter Paper I and Paper II), respectively. It was found 
that 6 Oph has a frequency spectrum which is similar to 
that of V836 Cen. In this paper we present our modelling 
based on accurate frequency determination and successful 
mode identitification obtained in Paper I and II. Our main 
objective is to test if the occurrence of core overshooting 
and non-rigid rotation found for V836 Cen also applies to 
e Oph. 
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The paper is organized as follows. In Section 2 we sum- 
marize the observational pulsation constraints which consti- 
tute the starting point of our seismic modelling of 6 Oph. In 
Section 3 we perform a detailed abundance analysis of 9 Oph 
with the aim to use the deduced metallicity as an additional 
constraint. In Section 4 we present the evolution and oscil- 
lation codes that we used in our study, as well as the physi- 
cal inputs. In Sections we derive the seismic constraints on 
global stellar parameters and on core overshooting. In Sec- 
tion 6 we test the hypothesis of a non-rigid rotation model 
thanks to the two observed multiplets. We end with a con- 
clusion in Section 7. 



Table 1. The pulsation modes derived from tlie photometric and 
spectroscopic results presented in Paper I and Paper II, respec- 
tively. The amplitudes of the modes are given for the u filter and 
for the radial velocities. 



ID Frequency (d ^) {i,Tn) u ampl. RV ampl. 

(mmag) (km s^-^) 





7.1160 


(2,-1) 


12.7 


2.54 


^2 


7.2881 


(2,1) 


2.1 






7.3697 


(2,2) 


3.6 






7.4677 


(0,0) 


4.7 


2.08 


I'S 


7.7659 




3.4 






7.8742 


(1,0) 


2.3 




1/7 


7.9734 


(1,1) 


2.4 





2 OBSERVATIONAL CONSTRAINTS 



The observational pulsation characteristics of 6 Oph derived 
in Paper I and II can be summarized as follows. The pho- 
tometric data (Paper I) were gathered in the framework of 
a three-site campaign, allowing the detection of seven pul- 
sation frequencies. The identification of their corresponding 
^- value showed the presence of one radial mode, one rota- 
tionally split £ = 1 triplet and three components of a rota- 
tionally split £ = 2 quintuplet. In addition, the spectroscopic 
observations (Paper II) lifted the ambiguity for the m- value 
of the observed £ = 2 main mode. The pulsation frequen- 
cies and their {£, m)-values are listed in Table[l] We note 
that such a frequency spectrum was observed for the star 



V836 Cen (Aerts et al. 2004 1. We also point out that the 
mode identifications of all the observed components of the 
quintuplet are determined for d Oph, which was not the case 
for V836 Cen. 

The position of appropriate models will be compared 
with the position of the star in the HR diagram that was 
determined photometrically and spectroscopically in Paper I 
and II. The obtained error boxes are represented in Fig.[l] 
We note that such a deviation between photometrically and 
spectroscopically derived effective temperatures is common 
for B-type stars (e.g. De Ridder et al. |2004[ Morel et al. 



20061 



Recently, Niemczura & Daszynska-Daszkiewicz ( 2005 1 
determined [M/H] for 6 Oph. However, they did not have 
information that 6 Oph is a triple system composed of a B2 
primary, a spectroscopic secondary with a mass lower than 
1 Mq (Briquet et al. 2005 I and a speckle B5 star (McAlister 
et al. |1993[ ). In what follows, we present a careful abundance 
analysis for the primary by taking into account the presence 
of the tertiary, the contribution to the lines of the secondary 
being negligible. 



3 ABUNDANCE ANALYSIS 

The non-local thermodynamic equilibrium (NLTE) abun- 
dances of He, C, N, O, Mg, Al, Si, S and Fe were calcu- 
lated using the latest versions of the line formation codes 
DETAIL/SURFACE and plane-par allel, fu lly line-blanketed 
Kurucz atmospheric models (Kurucz 1993 1. Curve-of-growth 
techniques were used to determine the abundances using the 
equivalent widths of a set of unblended lines measured in a 
mean CORALIE spectrum (see Paper II), which was created 




4.38 4.36 

log Teff 

Figure 1. The error boxes represent the position of 8 Oph in the 
HR diagram as derived from photometric (full line, Paper I) and 
spectroscopic data (dashed line. Paper II; dotted line, this paper). 
The positions of the models which fit exactly the three indepen- 
dent modes are also shown for the different couples (Z, «ov) given 
in TableH 



by co-adding the 86 individual exposures (all put in the lab- 
oratory rest frame prior to this operation). The reader is 
referred to Morel et al. ( 2006 1 for complete details on the 



methodology used to derive the elemental abundances. 

To correct for the contamination of the spectrum by the 
tertiary, our study is based on synthetic, composite spec- 
tra assuming the following parameters for this companion: 
Tcff = 19 000 K, log (7=4.0 dex [cgs] (Paper II) and a microtur- 
bulent velocity, ^=5 km s~^, typical of B-type dwarfs. We 
also considered that the tertiary contributes to 22% of the 
total light of the system in the optical band (Paper I). In ad- 
dition, we assumed the star to have a chemical composition 
typical of OB dwarfs in the solar vicinity (Daflon & Cunha 
2004). For iron, we assumed an abundance log e(Fe)=7.3 dex 
(Morel et al. 20061. We will examine below the sensitivity 



of our results to these assumptions. 

A standard, iterative scheme is first used to self- 
consistently derive the atmospheric parameters: Teff is de- 
termined from the Sill/III ionization balance, logg from 
fitting the collisionally-broadened wings of the Balmer lines 
and ^ from requiring the abundances yielded by the O II 
features to be independent of the line strength. We ob- 
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Table 2. Mean NLTE abundances (on the scale in which 
log€[H]=12) and total l-tr uncertainties. The number of used 
spectral lines is given in brackets. For comparison purposes, we 
provide in the last two columns the standard solar composition 
of Grevesse & Sauval (1998; Sun 1-D) and updated values in the 
present day solar photosphere derived from three-dimensional hy- 
drodynamical models (Asplund et al. 2005; Sun 3-D). 





e Oph 


Sun 1-D 


Sun 3-D 


He/H 


0.066±0.026 (10) 


0.085±0.001 


0.085±0.002 


log€(C) 


8.32±0.09 (7) 


8.52±0.06 


8.39±0.05 


log£(N) 


7.78±0.10 (23) 


7.92±0.06 


7.78±0.06 


loge(O) 


8.58±0.26 (27) 


8.83±0.06 


8.66±0.05 


loge(Mg) 


7.49±0.15 (2) 


7.58±0.05 


7.53±0.09 


loge(Al) 


6.24±0.14 (4) 


6.47±0.07 


6.37±0.06 


log€(Si) 


7.04±0.22 (8) 


7.55±0.05 


7.51±0.04 


log€(S) 


7.22±0.27 (5) 


7.33±0.11 


7.14±0.05 


log e(Fe) 


7.41±0.17 (27) 


7.50±0.05 


7.45±0.05 


Z 


0.0114±0.0028 


0.0172±0.0012 


0.0124±0.0007 



tain: Tcff =25 OOOilOOO K, log5=4.10±0.15 dex [cgs] and 
^=4j;3 km s~^. For comparison, we obtained in Paper II, 
Teff=24 000±1000 K and logff=4.1±0 .1 dex [cgs] using the 
NLTE code FASTWIND (Puis et al. 



20051. Other studies 



Dworetsky 1995 Morel et al. 2006 Kaiser 



show that using different methods on the same dataset can 
indeed lead to uncertainties of order 500 K for B-type stars 
(e.g. Smalley 
[20M|. 

The abundances are given in Table [2] and are com- 
pared with the standard solar mixture of Grevesse & Sauval 



(19981 and with values derived from time-dependent, three- 



dimensional hydrodynamical models (Asplund et al. |2005[ 
and references therein). The quoted uncertainties take into 
account both the line-to-line scatter and the errors arising 
from the uncertainties on the atmospheric parameters. Note 
that a possible downward revision of Toff by ~1000 K (see 
above) is explicitly taken into account in the total error bud- 
get. We infer a low helium content, but this quantity is un- 
certain and may be considered solar within the large error 
bars. On the other hand, there is no indication for the nitro- 
gen excess occasionally observed in other slowly-rotating (3 
Cephei stars (Morel et al. [2006). The resulting metallicity, 
Z=0.0114±0.0028, is identical, within the errors, to the most 
recent and likely realistic estimates for the Sun (Table [2|. 
To calculate this quantity, the abundances of the elements 
not under study were taken from Grevesse & Sauval ( 1998[ ). 
Our assumed neon abundance is indistinguishable from re- 
cent values derived for a sample of B-type stars in the Orion 
association (Cunha et al. 20061. The other species are trace 



elements and contribute only negligibly to metal mass frac- 
tion. 

To examine the sensitivity of our results to the vari- 
ous assumptions made about the physical properties of the 
companion, we have repeated the abundance analysis after 
varying the adopted effective temperature, surface gravity, 
chemical composition, and luminosity of the tertiary within 
the range of plausible values. Namely, we assumed in turn: 
Tcft—21 000 K, logg=3.7 dex [cgs], the abundances of all the 
metals enhanced by 0.3 dex relative to solar and a contri- 
bution of only 18% to the total light of the system in the 
optical, while keeping the other parameters unchanged. As 



Table 3. Sensitivity of the derived metal abundances and metal- 
licity of 6 Oph on the assumed properties of the tertiary. We 
quote the abundance differences compared with the values listed 
in Table [H 







-1-2000 K 


Alog g= 
-0.3 dex 


A log e= 
+0.3 dex 


flux ratio= 
18% 


A log 


.(C) 


-0.05 


-0.02 


-0.05 


-0.01 


A log 


e(N) 


-0.05 


-0.02 


-0.03 


-0.02 


A log 


e{0) 


-0.06 


-0.02 


-0.03 


-0.03 


A log 


e(Mg) 


+0.06 


+0.02 


-0.15 


+0.02 


A log 


e(Al) 


-0.05 


-0.01 


-0.04 


-0.01 


A log 


eiSi) 


+0.05 


+0.01 


-0.15 


+0.03 


A log 


e{S) 


+0.12 


+0.00 


-0.44 


+0.05 


A log 


e(Fe) 


-0.04 


-0.02 


-0.04 


-0.02 


AZ 




-0.0007 


-0.0003 


-0.0011 


-0.0003 



expected, the abundances of the chemical elements deter- 
mined from lines of low-ionization ionic species (e.g. Mg, S) 
are most strongly affected by the choice of the parameters 
for the cool component (Table [sf. However, the metallicity 
remains largely unaltered in all cases. Our conclusions re- 
garding the metal content of 6 Oph appear therefore robust 
against the exact nature of its speckle companion. 



4 STELLAR MODELS 

The numerical tools and physical inputs used in our study 
are the following. The stellar models for non-rotating stars 
were computed with the evolutionary code CLES ( Code 



Liegeois d'Evolution Stellaire, Scufiaire et al. 2007a I. We 
used the OPAL2001 equation of state (Rogers & Nayfonov 
|2b02| and Caughlan & Fowler (jl988| nuclear reaction rates 
with Formicola et al. (|2004[) for the "N(p,7)^^0 cross- 



section. Convective transport is treated by using the classical 



Mixing Length Theory of convection (Bohm-Vitense 19581. 
As shown in the previous section, the abundances of 9 Oph 
are in full agreement with the solar values of Asplund et 
al. ( 2005[ ). For the chemical composition, we consequently 
used the solar mixture from these authors, except for Ne. 
For this latter element, a direct abundance determination 
in a small sample of nearby B stars using photospheric lines 
(Cunha et al. |2006j ) suggests a value ~0.3 dex larger than 
quoted by Asplund et al. (20051. For our computations, we 



consequently adopted log e(Ne)=8.11. We used OP opacity 
tables (Seaton 2005 ) computed for the mixture in Cunha et 



al. ( 2006 1 (that is the mixture of Asplund et al. ( 2005 1 and 
loge(Ne)=8.11). These tables are completed at logT < 4.1 
with the low temperature tables of Ferguson et al. ( 2005 1 for 



the Asplund et al. ( 2005 1 mixture, as the effect of increasing 
Ne on low temperature opacities can be neglected for such 
a hot star. We computed stellar models with and without 
taking into account microscopic diffusion. For models with 
diffusion, we used the formulation of Thoul et al. (19941. 

Stellar models are parametrized by the initial hydrogen 
abundance X, the core convective overshooting parameter 
Qov, the metallicity Z, the mass M and the central hydrogen 
abundance Xc (which is related to the age). 

For each stellar model, we calculated the theoretical fre- 
quency spectrum of low-order p- and g-modes with a degree 
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of the oscillation up to ^ = 2 using a standard adiabatic 



code (Scuflaire et al. 2007b I, which is much faster than a 



non-adiabatic code but leads to the same theoretical pul- 
sation frequencies within the adopted precision of the fit, 
which was 10"'^ d~^. Once the models fitting the observed 
modes are selected, we checked the excitation of the pulsa- 
tion modes with the linear non-adiabatic code MAD devel- 
oped by Dupret et al. ('2001). 

In an attempt to explain the asymmetries of the ob- 
served multiplets, we also computed the adiabatic frequen- 
cies with the code FILOU (Tran Minh & Leon |1995| l, which 
includes effects of rotation up to the second order, according 
to the formalism of Soufi et al. ( |1998[ ). 



5 CONSTRAINTS ON STELLAR 
PARAMETERS AND CORE 
OVERSHOOTING 

5.1 Effects of diffusion 

9 Oph is a slow rotator with an equatorial rotation velocity 
of about 30 km s~^ (Paper II). Moreover, its surface convec- 
tion zone is very thin. In such conditions, diffusion mecha- 
nisms can occur and alter the photospheric abundances. We 
consequently investigated if diffusion could be the explana- 
tion for the marginal lower He content of 8 Oph compared 
to the solar value. Moreover, we checked its effect on the 
oscillation frequencies. 

Our calculations include microscopic diffusion (without 
radiative forces and wind, and using TBL94's routine; Thoul 



et al. 19941 and a turbulent mixing consistent with the re- 
sults of Talon et al. (19971. We were able to reproduce the 



observed surface metallicity and helium abundances with 
models having initially the solar composition {X — 0.7211, 
Y = 0.264, Z = 0.01485). Those models are very close to 
those obtained without diffusion, because the diffusion only 
affects the very superficial layers of the star. In particular, 
the models calculated with and without diffusion have ex- 
actly the same frequency spectrum. In Fig. [2] we show the 
metallicity and the Y/X profiles for a model which fits the 
three observed frequencies and has both a solar initial com- 
position and a surface composition compatible with the ob- 
servations. The diffusion and the turbulent mixing only af- 
fect layers down to a radius of 0.92 R. 

Even though these results agree very well with the ob- 
servations, it might be hazardous to trust them blindly. In- 
deed, radiative accelerations and stellar winds are very im- 
portant in those stars and can strongly affect the surface 



abundances (Bourge et al. 2006 1. They were ignored here be- 



cause calculations involving those effects are computation- 
ally intensive and difficult to perform, and clearly beyond 
the scope of this paper. 



5.2 Seismic analysis 

Since taking into account diffusion or not does not affect the 
derived stellar parameters of our models, we continued our 
analysis without diffusion but considering sufficiently large 
ranges for X and Z. 

For our seismic analysis we first searched for models 
that fit the radial mode with frequency 1/4 together with the 
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Figure 2. Metallicity and Y/X profiles for a model which fits the 
three observed frequencies and has both a solar initial composi- 
tion and a surface composition compatible with the observations. 
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Figure 3. The M — Z relations obtained by matching the radial 
mode and the central peak of the triplet for several values of the 
core overshooting parameter, X being fixed to 0.71 (dashed lines) 
or 0.7211 (full lines). 



zonal £ — 1 mode with frequency j/g. We then made use of 
the quintuplet to add additional constraints. 

We found that the radial mode is either the fundamental 
mode or the first overtone. However, the models with 1/4 as 
the first overtone are further away from the observational 
position in the HR diagram than the models with ^4, as 
fundamental. Moreover, none of the three modes is excited 
by the classical k mechanism for models with 1/4 as the first 
overtone, even for a value of Z of 0.015 with Y of 0.0264. 
We consequently concluded that the radial mode is identified 
as fundamental and a scan of stellar parameter space also 
reveals that the triplet is identified as pi. 
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Table 4. Relation between Z and Oo 
of the three independent modes. 



imposed by the matching 



z 


Ctov 


0.009 


0.51 


0.010 


0.48 


0.011 


0.45 


0.012 


0.43 


0.013 


0.40 


0.014 


0.38 


0.015 


0.37 



Table 5. Physical parameters of the model that matches the 
observed modes, X e [0.71,0.7211] and Z e [0.009,0.015]. 



M (Mq) 
Teff (K) 
logs (dex) 



8.2±0.3 
22260±280 
3.950±0.006 
0.38±0.02 
0.44±0.07 



Fitting one frequency results in finding one model 
along the evolutionary track for every combination of 
{X, Qov, Z, M). Then, X being fixed, fitting two frequen- 
cies gives a relation between two parameters for given values 
of the last one. For the considered values of cvov one thus gets 
several M — Z relations that are shown in Fig.jS] We can see 
that the fitting of the two frequencies implies an increase 
in mass if either the metallicity or the core convective over- 
shooting parameter increases. This figure also illustrates the 
order of magnitude difference in mass induced by different 
adopted X values. 

Because the zonal mode of the quintuplet was not ob- 
served, we computed it as the average frequency of the two 
surrounding modes with frequency and 1^2- A scan of stel- 
lar parameter space shows that the quintuplet is identified 
as g\. The fitting of a third independent frequency implies 
a relation between the metallicity and the core overshooting 
parameter: the lower the metallicity the higher the over- 
shooting. This Qov — Z relation is presented in Table|4] 

The positions in the HR diagram of the models repro- 
ducing the observed modes are shown in Fig.[l] for several 
couples [Z, Oov). One can see that the derived models are sit- 
uated in the cooler part of the photometric observed error 
box, and outside the spectroscopic one. This is actually the 
case for all studied j3 Cephei stars up to now and needs to be 
futher investigated. By considering a wide range of metal- 
licities Z G [0.009,0.015], for X G [0.71,0.7211], one obtains 
a core overshooting parameter Qov G [0.37,0.51] and a mass 
M G [7.9, 8.5] Mq. The other physical parameters are given 
in Table[5] The three modes are well excited by the classical 
K mechanism for a metallicity larger than 0.011. However, 
radiatives forces on iron allow this element to accumulate 
in the excitation region and lead to the excitation of addi- 
tional modes in lower metallicity stars (Pamyatnykh et al. 
2004[ Bourge et al. |2006[ ). Note that the only other theoret- 
ically excited mode with ^ 2 is the /-mode with £ = 1 but 
only for Z larger than 0.013. 



We point out that models matching the observed modes 
but computed using the solar abundances of Grevesse & 



Sauval ( 1998 1 with OPAL opacity tables are not excited for 
Z ~ 0.01. We refer to Miglio, Montalban & Dupret (|2006| 
for a detailed discussion on the implication of the adopted 
opacity tables and metal mixtures on the excitation of pul- 
sation modes. 

We also mention that the amount of overshooting found 
for 6 Oph corresponds exactly to that computed by De- 



upree (20001 by means of 2D hydrodynamic simulations of 
zero-age main-sequence convective cores. Our derived value 
is also in agreement with results obtained by Ribas et al. 



(20001 who provided an empirical calibration of convective 
core overshooting for a range of stellar masses by studying 
eight detached double-lined eclipsing binaries. They found a 
systematic increase of the amount of convective overshoot- 
ing with the stellar mass, the values being 0.3-0.6 for ~10 
Mq stars. 



6 CONSTRAINTS ON THE ROTATION 
6.1 First order analysis 

W hen the rotation frequency is small compared to both 
R^/GAI and the considered pulsation frequency, the pul- 
sation frequencies of modes differing only by the m value 
of the spherical functions are linked through a simple rela- 
tion. If we cissume that the rotational frequency Urot is a 
function of the radius r only, this relation reads 



= 1^0 +m I K{x)urot{x) dx , 
'o 



(1) 



where x — r/R. The rotational kernel K{x) depends on the 
considered mode (see Lynden-Bell & Ostriker 1967 or Unno 
et al. p89l ). 

The i = 1, pi triplet and three components (correspond- 
ing to m = —1, 1 and 2) of the £ = 2, gi quintuplet are ob- 
served. The components of these multiplets are not strictly 
equidistant as required by equation ([TJ. These departures 
from equidistance may result from the fact that the rotation 
velocity of the star is too large for a linear approximation 
to be valid. It may also result from the fact that a mag- 
netic field contributes to the splitting. As the departures 
from equidistance are not too large, we tentatively interpret 
them as errors in the measure of the splitting. So, we have 
Avi = 0.10375±0.005 d"^ for mode £ = 1, pi and Av2 = 
0.08457±0.003 d"^ for mode £ = 2, gi, denoting by Aiy the 
integral term in equation ([T|. 

Fig. [4] shows the behaviour of the rotational kernels for 
the modes we are investigating. It is clear that they do not 
probe the convective core. This was already the case for 



V836 Cen (Dupret et al. 2004 1 and, as in this case, with just 



two pieces of information on the behaviour of the rotation 
velocity, we analyse its trend inside the envelope by fitting 
the linear expression 



I^rot{x) = Urot,0 + Vrot,l{x ~ 1) . 

The splittings are then given by 

AVk = CkOVrotfi + CklVrot,! fc = 1, 2 



(2) 
(3) 
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Figure 4. The kernels for the £ = I, pi mode (soHd Hne) and for 
the £ = 2, gi mode (dashed line). The vertical dotted line marks 
the position of the boundary of the convective core. 



Table 7. Theoretical rotational splittings (d"'^) obtained with 
a second order perturbative treatment of rotation, for a model 
fitting the observed zonal modes. 



Solid rotation: u^ot = 0.10751 d ^ 

^=1 1^1,0 - I'l.-l i^i.i - ^1,0 

0.10934 0.09816 

^ = 2 (1^2,1 - i^2,-l)/2 1^2,2-1^2,1 

0.08661 0.08859 

Sohd rotation: u^-ot = 0.10682 d-^ 

^ = 1 1^1,0 - fl.-l 1^1,1 - 1^1,0 

0.10860 0.09756 

^ = 2 (1^2,1 - i^2,-l)/2 !^2,2 - 1^2,1 

0.08605 0.08800 

Differential rotation: frot = 
^=1 1^1,0-1^1,-1 

0.10932 
i = 2 {V2,i-v2,-i)l2 

0.08605 



1^1,1 - Vl-l 

0.2075 

!^2,l - !^2,-l 

0.1732 



1^1,1 - Vl,-1 

0.2062 

1^2,1 - i^2,-l 

0.1721 



0.10915 + 0.00549 {x - 1) d" 
1^1,1 - 1^1,0 1^1,1 - i'i,-i 



0.09819 

^^2,2 — !^2,1 

0.08826 



0.2075 

!^2,l - i^2,-l 

0.1721 



Table 6. Observed rotational splittings (d ^). The subscripts of 
vo m denote the degree £ and azimuthal order m. 



^=1 I'l.O - 


1^1,1 - i^l.O 


1/1,1 - 2^1,-1 


0.1083 


0.0992 


0.2075 


£ = 2 {y2,i-V2,-i)l2 


!^2,2 — 1^2,1 


1^2,1 - 1^2,-1 


0.08605 


0.0816 


0.1721 



with k referring to the two known splittings, and 



CkO 



Ckl 



Kk (x) dx , 



Kk{x)(x — l)dx . 



(4) 
(5) 



Taking the errors on the splittings into account, the sys- 
tem ([3| gives (all frequencies in d~^) 
0.0878 <Urot,o < 0.1419, 
-0.0491 < Vrot,i < 0.0971 . 
The splitting data are thus consistent with a constant rota- 
tion velocity inside the envelope, a rotation period of 9.2±2.2 
days and an equatorial rotation velocity of 29±7 km s~^. 
This latter value is in full agreement with the v sin i and 
equatorial rotation velocity derived in Paper II. 



6.2 Second order analysis 

A small asymmetry is observed in the £ = 1 and 1 = 2 
multiplets, as shown in Table [6] It is well known that asym- 
metries are explained by the effect of terms of higher order 
in i-'rot /j/puis in the pulsation equations. We determine here 
the adiabatic frequencies with the code FILOU (Tran Minh 
Suarez 2002 1 . In the version of the code used 



& Leon pM 

here, the effects of rotation are included up to the second 
order, following Soufi et al. ( |1998[ ). This code needs as in- 
put the spherically symmetric component of the structure 
model; it determines a posteriori the second order defor- 
mation due to rotation. In principle, the gravity must be 
corrected for the effect of centrifugal acceleration already 



in the spherically symmetric component of the model. This 
correction would affect very slightly the frequencies (slow 
rotation) and has a negligible effect on the multiplet asym- 
metries. It is not included here, which allows us to use as 
input one of our best no-rotation model fitting the zonal 
mode frequencies, as determined in the previous sections. 
The main global parameters of this model are: M = 8.4 Mq, 
log(L/L0) = 3.7346, Tcs = 22053 K, X = 0.72, Z = 0.014 
and aov = 0.38. 

In Table [7] we give the theoretical rotational splittings 
obtained with this model and the second order treatment of 
rotation. We recall that the second order terms cancel in the 
combination i'e,m — i'i,-m- The last column gives such com- 
bination, which we use in the fitting procedure. Comparing 
columns 1 and 2 shows the splitting asymmetry. For the first 
results given in this table, we consider a rigid rotation. In 
the first case, the rotation frequency is u^^t ~ 0.10751 d~^. 
With this value, we fit exactly the observed value of the 
£ = 1 splitting: i^i^i — vi.-i = 0.2075 d~^. In the second 
case, the rotation frequency is v^-ot = 0.10682d~^. With this 
value, we fit exactly the observed value of the £ = 2 splitting: 
1^2,1 — 1^2,-1 = 0.1721 d~^. In the last case of Table [t] we 
consider a differential rotation law of the same linear form 
as Eq. ([2J. The coefficients are adjusted to fit at the same 
time the observed £ = 1 and 1 = 2 splittings 1^1,1 — 1^1,-1 
and 1/2,1 — i^2,-i. This gives the linear differential rotation 
law: Urot = 0.10915 + 0.00549 (a: - l)d-^ As for the first 
order analysis, we see that rigid rotation models cannot be 
eliminated. 

The results of Table [7] show that the non-equidistance 
of the £ = 1 triplet is relatively well reproduced by the sec- 
ond order effect of rotation. However, the asymmetries of 
the £ = 2 multiplet do not fit the observations. This discrep- 
ancy could come from observations (multiplet not entirely 
resolved) or theory (effect of higher order terms). We did 
the same analysis with other models fitting the zonal mode 
frequencies and find quasi-identical results for the rotation 
velocity. 
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Our study of the /3 Cephei star d Ophiuchi is a new illus- 
tration of the power of asteroseismology for this class of 
pulsators. The couple {X, Z) being chosen, the observation 
of three independent modes is enough to derive the other 
parameters that characterize a stellar model, for adopted 
physical inputs. We point out that such a success can also be 
attributed to the unique derivation of {£, m) thanks to state- 
of-the-art empirical mode identification techniques used in 
Paper I and II. 

A detailed NLTE abundance analysis showed that the 
considered abundance values and thus the metallicity of 
6 Oph correspond to the new solar mixture of Asplund et al. 
(20051. In particular, the CNO abundances are much more 



consistent with the 3D values of Asplund et al. ( 2005 I than 
with the ID values of Grevesse & Sauval (19981. This is 



generally the case for B-type stars (Morel et al. 20061. We 
found a mass M = 8.2±0.3 Mq and a central hydrogen 
abundance Xc = 0.38±0.02 for the star. 6 Oph is the fifth 
P Cephei star for which the occurrence of core overshooting 
is deduced by seismic interpretation and is the target with 
the highest derived value {uov = 0.44±0.07) among them. 
However, it might be that the core overshooting parameter 
of previously modelled (3 Cephei stars is underestimated, 
as the case of V836 Cen illustrates it. For the modelling of 
this star, Dupret et al. (20041 adopted a value of Z larger 



than 0.016 in order to get the excitation of the modes. How- 
ever, Morel et al. (2006,) determined Z = 0.0105±0.00 22 for 
V836 Cen. In addition, Miglio, Montalban & Dupret ( |2006 1 
recently showed that modes can be excited for Z ~ 0.01 if 
one uses the new solar abundances together with the OP 
opacities, which was not the case in Dupret et al. (2004 1. 



Finally, considering a lower value of Z increases the Oo 
the star (see Fig. 3 in Dupret et al. 



of 



20041 



We also showed that the asymmetry observed in the 
£ = 1 triplet can be well reproduced by taking into account 
the effects of rotation up to the second order. For the quin- 
tuplet, the agreement is however not as good. Contrary to 



V836 Cen (Dupret et al. 2004 1 and u Eri (Pamyatnykh et al. 
|2004 l) for which non-rigid rotation was proven, the observed 
rotational splitting of two modes for 9 Oph are still com- 
patible with a rigid rotation model. In the near future, we 
can expect stronger constraints on the internal rotation of 
/3 Cephei stars from data collected from space missions (e.g. 
MOST, COROT). With the observation of rotational split- 
ting of many modes having different probing kernels, we aim 
to determine the internal rotation law of massive B-type 
stars. 
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